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Abstract 

We derive non-perturbative sum rules in SU(A^) lattice gauge theory at finite tempera- 
ture. They relate the susceptibilities of the trace anomaly and energy-momentum tensor 
to temperature derivatives of the thermodynamic potentials. Two of them have been de- 
rived previously in the continuum and one is new. In all cases, at finite latttice spacing 
there are important corrections to the continuum sum rules that are only suppressed by the 
bare coupling g^. We also show how the discretization errors affecting the thermodynamic 
potentials can be controlled by computing these susceptibilities. 



1 Introduction 



Sum rules in continuum QCD at zero temperature were introduced by Novikov et al. [T] 
and a lot of hadron phenomenology was subsequently based on them. A few years later 
Michael derived sum rules for SU(iV) pure gauge theories in lattice regularization [2l[3lll], 
which were only recently generalized to Wilson lattice QCD in [^j. Since sum rules are 
relations that hold non-perturbatively, the lattice regularization provides a framework in 
which their derivation proceeds in a particularly rigorous way: it only involves operations 
on multi-dimensional integrals. 

On the lattice the simplest identities relate zero-momentum three-point functions to 
the spectrum of the theory. By comparing the sum rules to continuum relations [6], one 
realizes [71 [8] that they relate the normalization of a particular discretization of the trace 
anomaly and the energy-momentum tensor to anisotropy coefficients. The latter are deriva- 
tives of the bare lattice parameters with respect to physical parameters such as the lattice 
spacing in hadronic units or the ratio of spatial and temporal lattice spacings. Indeed this 
normalization is non-trivial since translation invariance is broken down to a discrete group 
at finite lattice spacing. 

Ellis et al. derived finite-temperature sum rules in pure gauge theories [9] and in full 
QCD [TU]. In this paper we rederive the SU(A^) gauge theory sum rules, focusing on those 
concerning two-point functions of the trace anomaly and the energy-momentum tensor, 
in lattice regularization. We find that they have important corrections to the continuum 
versions, which are only suppressed by one power of the bare coupling g^. From the point 
of view of Monte-Carlo simulations, where thermodynamics calculations are performed 
around ^ 1, they can thus not be neglected. We also derive a new sum rule involving 
only the traceless part of the energy-momentum tensor, and we relate the results obtained 
to contact terms in the two-point functions of the Hamiltonian. 

The p = two-point function of the trace anomaly, in other words its susceptibility, is 
related to the rate of change of (e — 3P) /T^ with temperature [S] (e is the energy density 
and P the pressure). Because the bulk viscosity is related to this two-point function by 
a Kubo formula [llj, it was argued recently [12] that the bulk viscosity rises sharply just 
above the deconfining temperature Tc. Direct calculations of the two-point functions at 
p = and general lo = have confirmed the existence of this effect |13j . In the context of 
such calculations, the sum rule can be used to constrain the reconstruction of the spectral 
function p{u}). 

Another application of these considerations to finite-temperature Monte Carlo simula- 
tions is to compute directly the leading lattice spacing dependence of the thermodynamic 
potentials. This idea has the most potential of being useful in the context of full QCD 
simulations, where the computational cost is high and grows with a large power of the 
inverse lattice spacing. 

Decomposing the energy-momentum tensor T^^, into a traceless part O^y and a scalar 
part 6 via T^^ = 6^^ + \5^i,9, the explicit Euclidean expressions are 

9{x) = f3{g)/{2g) F-^{x)F-^{x) e,,{x) = \5^,F-„F-^ - F^^F^^. (1) 

The beta-function is defined by qdg/dq = P{g) = -f{bo+bif+. . .) and 6o = llA^/(3(47r)2), 
bi = 34A^^/(3(47r)^) in the SU(A^) pure gauge theory. The gauge action reads \F^^F^j^ in 
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this notation. If (. . .)t denotes the thermal average at temperature T, 

e-3P={e)T-{e)o, e + P = l{9oo)T. (2) 

In section 2 we introduce our notation, review the relations relevant to thermodynamics 
and introduce new anisotropy coefficients. In section 3 we derive the sum rules on the 
lattice. In section 4 we take the extreme continuum limit, gfg <C 1 and compare our results 
to those of [9j . Section 5 describes the possibility of computing the leading discretization 
errors affecting e and P in numerical simulations, and section 6 contains concluding remarks. 



2 Thermodynamics and the energy-momentum tensor 

We consider a Euclidean lattice of spatial extent N^- sufficiently large that the thermody- 
namics limit has been reached. The time-extent A^,- = Lq/u fixes the temperature T = I/Lq. 
Thermal averages are denoted by (...). The temperature dependence is made explicit by 
(. . .)t to distinguish this average from the average (. . .)o on a zero-temperature lattice. 
The contents of this section is mostly a review of the early lattice thermodynamics arti- 
cles [14[ riB] , with the exception of subsection 12. 4i We however emphasize a lot more the 
role of ^00 and 9, since they are the operators of interest in the sum rules. 

2.1 Isotropic lattice 

We start from the Wilson action |18j for SU(A'^) gauge theories: 

-^g = /3EEV(^)> (3) 
S^^x) = ^ReTT{l-U^{x)U,{x + ail)U^{x + aiy)-'U,{x)-'} (4) 
and (3 = It is useful to consider two separate sets of parameters: 

bare parameters: /3, A',- (5) 

physical parameters: o.{(3), T{(3,Nr) (6) 

where 1/T = Lq = NT-a{0). With the notations 

S± = Sa±Sr, Sa = ''^Skl, 5t- = ^SoA:, (7) 

k<l k 

we use the following discretizations: 

G(x) = S+, eoo(x) = Z-{P) S., (8) 

with 

^'(« = 5i|^ (9) 

The presence of the normalization factor Z{P) = 1 + 0(1//?) must be expected, since 
f d^Ji.6oo{x) is not a Noether charge, due to the lack of continous translation invariance on 
the lattice. A precise expression can be given for Z{p) in terms of derivatives with respect 
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to the anisotropy, see Eq. [29]and also Eg. 1581 For a parametrization of Z{P) in the case of 
SU(3) based on the data of [17J, see [16j (Eq. 6). The continuum Umit takes the form 

F^{P,N,)^N^mT-{e)o) '-=^^U{T) (10) 

F_(/3,iV^)^^iV,4(eoo) -° '-^^f-iT), (11) 
where the leading corrections are O(a^). 



2.2 Anisotropic lattice 

On the anisotropic lattice with spatial lattice spacing Oo- and temporal lattice spacing Ur, 
the action reads 

Sg = ^PMx)+PrSr(.x). (12) 

X 

The two separate sets of parameters are: 

bare parameters: Pa, Pt, (13) 

physical parameters: aa{Pa, Pr), (.{Pa, Pr), T{P„, (3r, Nr) (14) 



where 



C = a^/ar (15) 
l/T = Lo = Nrar = Nra„C^. (16) 

Obviously, ^ = 1 when P^j = P^-. 

We use the following discretizations: 

r'G(x) = Z+{Pa,pr)Sa + Z+{Pa,Pr)Sr, (17) 
r'eoo(x) = Z-{P,,Pr)Sa-Z-{Pa,Pr)Sr, (18) 

where at the symmetric point ^ = 1, 

Z+(/3,/5) = Z+(/5,/5) = Z+(/3), Z-{P,P) = Z;iP,P) = Z'iP). (19) 

The continuum limit Oo- ^ is taken at fixed ^. The factor Z^^ are such that, for 
instance, 'S>oo{x))^{J d'^x6oo{x)). The continuum limit of thermodynamic potentials 
is obtained according to 

F+{Pa,Pr,N^)^N^rMQ)T-{e)o) '-^^U{T) (20) 

F.{Pa,Pr,N^)^^KcH^m) ^^/-(r), (21) 

where the leading corrections are O(a^). 
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2.3 Thermodynamics and normalization of 9 and O^j,,, 

In this section we relate the normahzation factors Z^^ to derivatives of the bare parameters 
with respect to physical parameters. We start from the thermodynamic relations 

1 d\ogz 1 d\ogz 

where 

log Z{P,, Pr,N„, Nr) = log Z(/3,, Pr,N,, Nr) " ^ log Nf). (23) 

The subtraction, which sets the free energy F = —T log Z to zero at a reference temperature 
^ref = / {^T^^ o,t) , is necessary in quantum field theory. On a ^ = 1 lattice a common choice 
is A^^'^^ = Nfj, which implies that T^ei = in the thermodynamic limit. We can combine 
the equations 

dlogZ dlogZ 

and -— = 



into 



dlogaa ^logC 



{e-3P)alar = -j^^{S„)t-o + ^^^{Srh-o (24) 
-^{e + P)a,ar = - 777— 7 + TTTT:-— ('5.) - hrr^ (25) 



. d log ^ 4 9 log Oo- / \d log ^ 4 9 log 

Prom here we read off the normalization factors of O and Oqo: 

= = (26) 

aioga^ ^ alogoo- 

^ aiogC 4aioga^' ^ ^ 51ogC 451oga^' ^ ' 

Since, by Euclidean symmetry, Z~ Z~ , we have the equalities 



dWa + Pr) ?=1 1 dp 



d log 2d log a 



(28) 



diPr-P.){a,,0 «=i 2;3Z(/3). (29) 
alog^ 

We discuss a different choice of bare parameters often used in numerical simulations in 
appendix A. 

2.4 Derivatives of Zj^ at ^ = 1 

At ^ = 1, {dp^ + dp^) becomes d/dp. Using 

( d log a<T d log gp- \ , / dPr dPa \ 

a/3,, a/3, 5=1 1 aiog^ aiog^ 

a log 5 aiogg ofl'7/'fl\ a/3r a/j^ > 

a/3a a/3r / '^I^^KI^JdT^ \ dlogaa dloga^ ) 



(30) 
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one easily obtains the relations 



{Zt + Zt) '= 3- 



2\dfia dPr'' ^ dloga/3Z(/3)' 

d_\fr,+ _^+. 5=1 dp dp{f3Z{f3)) 

2\dPa dl3r)^ " ^' dloga PZ{P) ' ^ ' 

We shall need these relations in the next section. Similarly we introduce the quantities 

,±,, , , 1/ a a 



A„V/s„A)-5(j^-a3;)K±2.-)- (32) 



At = 1 they evaluate to 



1 dp 



1 dZ 



2(ilogaL/3 Zdp 



1 d\p^-pr) 



= 3-o:7TZ7Z^ + ^ + (33) 



2/3Z(/3) d[\og£,f 



2L 8a(loga)"^ a(log^)^ J 

These derivatives thus depend on second derivatives with respect to ^. In appendix B we 
obtain the leading order values of A(]q(/3) in g^. 

3 Derivation of the sum rules 

We now derive the sum rules, neglecting O(a^) discretization errors, but without using 
perturbative approximations to normalization factors such as ^fj^^^ and Z{P). 

3.1 Derivation on the isotropic lattice 

We consider a renormalization group invariant (RGI) quantity f{a,T), which is obtained 
as the continuum limit of a function -F(/3, A'^^) of the bare parameters. The renormalization 
group equation adf/da = implies 

dT dloga dp ^ ' 

We have used 

adNr/da = -N^ and N^dN^F = -Tdrf. 

In particular, we can apply this equation to F±(P, Nr), since they are RGI quantities (see 
Eq. [lOl E]). For the case of F+, we obtain 

a-^E.0(^)©(O))^-a-^E.©(^)0(O))§ = T^5T^ + ^,37ilHi^ (e-3P). (36) 

This sum rule was first derived in [9J in the continuum, in which case the second term on the 
right-hand side is absent. Indeed, the factor multiplying (e — 3P) behaves asymptotically 
as 26ig(Q at small bare coupling. Consider next the case of the sum rule reads 

^ (E. e(x)eoo(0))^ = T^Ot'-^ + {e + n (37) 

Note that the left-hand side vanishes by Euclidean symmetry at T = 0. The factor multi- 
plying (e -|- P) in the second term on the right-hand side vanishes in the continuum limit 
as (-26o9o)- 



5 



3.2 Derivation on the anisotropic lattice 

For any RGI quantity f{aa,£,,T), Uada^f = and S.d^f = respectively imply 

T%(~l) = (^ (38) 

We have used 

aada^ Nr = -Nr and N^dN^FiP^ ,Pr,Nr) = -Tdrf. 

From now on we evaluate the expression at the isotropic point (3^ = (3^- The determinant 
of the matrix is then 

A = 2/3Z(/3)-^. (39) 
a log a 

Taking suitable linear combinations, we obtain the two equations 

_T|i = J!!-(?L^?e:\ (40) 

The first relation is equivalent to Eq. [33] derived on the isotropic lattice, since = 
{dy + dz)f\y^z=x for a general function of two variables {y,z). We therefore focus on the 
second relation in the following. 

The observables f± are RGI quantities. Consider first f+{T). Using Eq. [3l] and the 
thermodynamic relations TdxP = e + P and (e — 3P)/T'^ = Tdxip/T^), Eq. ST] leads to 
Eq. [37] derived on the isotropic lattice. 

We now apply Eq. [Uto f-{T). We obtain a new sum rule, 

a'(E.©oo(x)eoo(0))5^-/3Z(/3)Aoo(/3)a-'(5+)T = |Ao+o(/3)(e + P) + (f) V^t^. (42) 

The quantities A^q(/3) are defined in Eq. [32] Since the right-hand side of Eq. [l2l manifestly 
has a finite continuum limit, this equation implies that the short-distance quartic diver- 
gence of the integrated correlator is compensated by the quartic divergence (a~^) of the 
expectation value of the trace anomaly, 

{PS+)t = ^dA (1 + 0(<7o')) {dA ^N^- 1). (43) 



3.3 Contact terms in two-point functions of the Hamiltonian 

The (06), {OqoOqq) and {OqoO) correlators are related at vanishing spatial momentum because 
the Hamiltonian operator f d^xToo has simple correlation functions: 

{Jd^^ roo(xo, x) = T^dT{0)T + Ao{T) 5{xq) (44) 

for any local operator O. The delta function arises because the Hamiltonian operator 
applied on transfer-matrix eigenstates with energies at the cutoff scale does not yield the 
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expected matrix elements; for instance off-diagonal matrix elements are expected to appear 
in general. 

The sum rules (Eg. [36 | \37 \ [i2|l determine the contact terms Ag^^ and Ag: 



An 



^00 (go )^ (go) 
gldg^'^ /d\oga' 

1 <f9o' 



(^>T + 3 iA+(5o) + 



go dg^ 
16 dlog a 



Zdgl 



1 (e + P),(45) 



4d(loga)2 dgQ'^/dloga 



3 2 dg, 



-2 



1 (e-3P) + -5o' 
' 4 d log a 



\_gldz_ 



(46) 



The contact term of (TqoToo) is then given by At^q = Aggg + ^Ag. Note that the con- 
tact terms have a quartically divergent contribution, plus finite, temperature-dependent 
contributions. 



4 Sum rules in the continuum 

Taking the bare coupling -C 1 in Eq. ([SBj [571 1121) yields the following continuum sum 
rules: 

f — 3P 

(Kxe(x)e(O))^-(Kx^(x)0(O))S = T'dT-^ (47) 

{fd'xe{x)9oo{Om = It'Bt'-^ (48) 

(Kx^oo(x)^oo(O))^ + 2^(0)t = \Ko{^ + P) + {^^T^dT'-^. (49) 

The coefficients A^q are now to be taken at g^ = 0, where they are pure, finite numbers. 
We compute these numbers in appendix B, see Eq. ([67]) and ([7^) . The calculation of A^q 
suggests that the latter is independent of the regularization used. If true, this would mean 
that the regularization dependence cancels entirely between the two terms on the left-hand 
side of this equation. It would be useful to derive Eq. HU] in a different regularization to 
confirm this. 

The difference of relation (|49p between finite and zero-temperature gives 

(Kx^oo(x)0oo(O))^-(Kx^oo(x)^oo(O))§ + ^(6-3P) = f A+o(6 + P) + (|) 

This relation shows that even after subtraction of the quartic divergence, a temperature- 
dependent logarithmic divergence remains in the susceptibility of Oqq. 

5 Sum rules and cutoff effects on e and P 

The idea to remove the leading cutoff effects on physical quantities by using lattice sum 
rules was proposed in [5j. Here we show that it can be applied to thermodynamic potentials. 
Consider for instance (e — 3P)/T^. On a = 1 lattice, this quantity is obtained by taking 
the Nt — > oo limit of 

ip{Nr)=FMNr),Nr), (50) 
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where f3{Nr) is tuned so that {Nt-o) is constant and was defined in Eq. [101 Following 
the steps of section 3, we can evaluate 

_ dF+ 
dlogNr ~ d log Nr 



dlogNr 

Thus the cutoff effects can be evaluated in Monte-Carlo simulations at fixed /?. The first 
term is itself unambiguous only up to O(a^) if a symmetric difference scheme is used, and 
0(a) if not [5]. It requires performing a simulation at a second value of Nr- Thus in 
total three simulations are required (for instance with the number of points in the time 
direction set to N^, Nr + 1, and N^j for the zero-temperature subtractions). Choosing a 
different couple {l3{Nr), N!^) tuned to the same temperature requires four simulations in 
total and provides essentially the same information (unless N!^ is much larger than N^-, but 
in practice, typical values are A',- = 6 and N!^ = 8). If one follows both strategies, one can 
check how close '^{Nlj.) is from 

If N-j- is large enough, ip is in the regime where O(a^) effects dominate over higher order 
cutoffe effects and '^{N'r) will be numerically consistent with this expression. In general, this 
provides a way of testing whether is in this regime without having to perform simulations 
at N'^ > N'^. Since the cost of finite-temperature calculations grows with a high power of 
A'^,-, this information is very precious. 

6 Conclusion 

We have derived finite temperature sum rules, valid at finite lattice spacing up to O(a^) cor- 
rections. The main results are Eq. (j361 137 1 H2]) . and, for the reader interested in continuum 
results, Eq. (fi7Hi9]) . 

As an application of these considerations, we have proposed a way to check whether 
thermodynamics calculations are performed in the regime where the 0{a?) cutoff effects 
dominate over higher order cutoff effects, using only two values of A^^^- 

Further sum rules can be obtained for other RGI quantities. Equations ESI and (j4U|41|) 
can for instance be applied to renormalized Polyakov or Wilson loops in order to study 
thermal contributions to quark masses, and the static potential relevant to J/^ suppres- 
sion [20j. Finally the sum rules can be generalized to full QCD with commonly used quark 
actions. 

This work was supported in part by funds provided by the U.S. Department of Energy 
under cooperative research agreement DE-FG02-94ER40818. 



d/3 dF^ 
d\oga d(3 



dilogaY 



(51) 
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A A different choice of bare parameters 



Although the set of bare parameters {f3a,(5r) is most convenient to derive sum rules, in 
numerical practice, it is more convenient to parametrize these parameters as 



(52) 



In order to take the continuum limit at fixed anisotropy ^, the first task of the lattice 
practitioner is to establish the lines of constant ^ in the (/?, i^o) plane, so that can 
thereafter be viewed as a function of {P,^,)- Secondly the relation between /3 and must 
be worked out at fixed anisotropy ^. After this preparatory work, the set of variables used 
in practice is (/?, ^). 

The expression (|26p can thus be written as 



9 log tta 

Slog a^j 



ia/3(a,,0r-, /9 9^0 (AO 



^0 91oga^ L 
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d log a„ 



1 

1 + 



Co dp 



(53) 
(54) 



Similarly, using 

5/3(a,,0 



9/3(a, , e) a log a, (/?, 5/3(a, , 5 log Co (/3, ^ log a, (/3, ) 



51og^ 
we obtain 

d(ia{a„,i) 
31og^ 



d log d log ^ d log Oo- 

^ aioggo(/?,0 ^ 
^0 91og^ 

9/3(a<^,0 9 log aa(/3,?o) 



51og^ 



91og^o 



9/3r(Qa,0 

91og^ 



1 + 



a log a log ^0 
5 log ^0 (AO 



gloggo(AO 
dp 



(55) 



(56) 



^ _ a/3(a^,0 91oga^(/3,eo) /l ^ 91ogeo(A6 
Sloga^ 51og^o 

8/3n-,T(aa-,g) 



These expressions suggest how to determine g[ogg non-perturbatively. Since, by Eu- 

— ' ^^'^op^^ ^= 0, we have the equalities 



clidean symmetry, 



and 



1 dp 



Z{P) 



4 d log o 



(57) 
(58) 



B Leading-order computation of A 



± 

00 



In this appendix we calculate the coefficients A^q(/3) defined in Eq. 
notation = 2sin(p^/2), we define the dimensionless integrals 



-y 



Po 



Pi 



Using the standard 
(59) 
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(2vr)3 eoPl + T.,pI 



(60) 



The variable pq takes the values 2'Kk/NT- for < A; < Nr- One finds that /(j(l,oo) 
cxd) = 1/4 and 



5^0 



1 



-4 



-0.154933... 



(61) 



B.l Aqo 

Notice first that \qq{(3) can be rewritten 

A„-„(/3) = iz-(,3)(|--|-)log(|). (62) 

The matrix elements of Sqo on physical states are RGI quantities. On a = 1 lattice, 

(fi|0oo|f^) =0, (63) 

as a consequence of the Euclidean symmetry on an Nfj = Nj- = oo lattice. Therefore Eq. [63] 
must be satisfied also on a 7^ 1 lattice. This condition determines the ratio /Z^: 



Z- (Sr) 



Zt {Sa)o 



At leading order on an A'^^ x 00'^ lattice, 



^{Sr)T = ^[lr{eo,Nr) + lAeO,Nr)]=l-3Ue„N,) 



(T = (a,iV,)-i), 



which leads to 



dl 

Vo = 1 + 8 T^iCi = 1, 00) = 0.586844 . 



(64) 

(65) 
(66) 

(67) 



B.2 A 



00 



A second physics condition (in addition to Eq. is necessary in order to fix Z^ and Z^ 
separately and therefore to determine AJq. We impose the condition 



3 

^(eoo)T = |(e + P)/T^ 



(68) 



This leads to the expressions 

z- = 
z- = 



3 e + P 
4dAT4 

3e + P 
AdAT^ 



Pa 



{(3rSr/dA)o 

{/3„Sa/dA)o 



(69) 
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where 

Wriil) = {PrSr/dA)0 WaSJdAh " {PrSr/dAh {(3^SJdA)o- (70) 

Cutoff effect due to finite A',- can be removed by taking the hmit A',- — > oo. Expressions 
Eq. [69]in principle allow for a non-perturbative determination of Z~^, but at tree level we 
shall use the Stefan-Boltzmann expression vr^dyi/lS for the right-hand side of Eq. [HSl In 
that approximation we have 

Wrifo) = I [lACl Nr) - UCl'^)] ■ (71) 
For = 1, we know that both Z~ and Z~ are equal to (3 at leading order. Thus 

hm N^Wt{1) = ^. (72) 

AT,-— »oo ZU 

Because limj\^^^co Nt^t{Co) is a continuum limit, p can be replaced by p and one then 
finds that 

hm Wrifo) = fo N^Wt{1) = ^. (73) 

One then straightfowardly obtains 

A+o = 6. (74) 
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